Epitaxial-strain-induced multiferroicity in SrMnOs from first principles 



0^ 
O 
O 
(N 

CJ 

Q 



o 

> 
oo 
cn 

d 

O 



X 



Jun Hoc Le^ and Karin M. Rabe 

Department of Physics and Astronomy, Rutgers University, Piscataway, New Jersey 08854-8019, USA 

First-principles density-functional calculations reveal a large spin-phonon coupling in cubic 
SrMnOs, with ferromagnetic ordering producing a polar instability. Through combination of this 
coupling with the strain-polarization coupling characteristic of perovskites, the bulk antiferromag- 
netic paraelectric ground state of SrMnOs is shown to be driven to a previously unreported multi- 
ferroic ferroelectric-ferromagnetic state by increasing epitaxial strain, both tensile and compressive. 
This state has a computed polarization and estimated Curie temperature above 54 fiC/cm^ and 
92 K. Large mixed magnetic-electric-elastic responses are predicted in the vicinity of the phase 
boundaries. 
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Multifcrroic materials have been the subject of con- 
tinuing attention both for fundamental physics and for 
potential apphcations including transducers and infor- 
mation storage [l|, [3, S] ■ There is particular interest in 
the search for multiferroic materials with large polariza- 
tion {P > 1/^C/cm^) and magnetization that persists to 
high temperatures, as well as a strong coupling between 
the magnetism and the electric polarization. This cou- 
pling may be enhanced by proximity of the system to a 
morphotropic phase boundary. 

With epitaxial strain, it is possible to widen the search 
to include phases and phase boundaries that do not ap- 
pear in bulk. In many paraelectric perovskite oxides, 
epitaxial strain couples strongly to the lowest-frequency 
polar phonon; this is responsible for the phenomenon of 
epitaxial-strain-induced ferroelectricity, which has been 
intensively studied both experimentally and theoretically 
[1, 0, Q. In an antiferromagnetic-paraelectric system, 
there is a further intriguing possibility. If the system 
has a spin-phonon coupling in which the lowest-frequency 
polar phonon is softer for ferromagnetic ordering than 
for antiferromagnetic ordering, then epitaxial strain en- 
hancement of a polar instability lead to the lowering 
of the energy of the ferromagnetic (FM)-ferroelectric 
(FE) state below that of the antiferromagnetic (AFM)- 
paraelectric (PE) state. This mechanism for producing 
a multiferroic phase was proposed and elucidated using 
first-principles calculations for EuTiOa , with polariza- 
tion above 30 /xC/cm^) at -2.0 % strain but a low Curie 
temperature comparable to the bulk Neel temperature of 
5.5 K [|. 

The demonstration of the spin-phonon coupling mech- 
anism for epitaxial-strain-induced multiferroicity in 
EuTiOa suggests a search for other paraelectric antifer- 
romagnetic materials in which this mechanism could be 
realized and the performance optimized. The primary 
criterion is that of a large downward shift in the frequency 
of the lowest polar phonon with ferromagnetic ordering, 
which leads to a large energy gain for polar distortion of 
the ferromagnetic state. Next, we look for systems with 
a moderate Neel temperature (T^)' while it woidd be 



preferable to have the magnetic ordering temperatures 
as high as possible, this is limited by the fact that the 
FM-AFM energy splitting cannot be larger than the scale 
of the energy gain for polar distortion of the ferromag- 
netic state. Finally, the polar instability should be strong 
enough to compete with any other structural distortions, 
such as oxygen octahedron rotations. 

i?-site magnetic perovskite oxides are of particular in- 
terest, as the larger exchange coupling results in much 
higher magnetic ordering temperatures than that of A- 
site rare-earth systems such as EuTiOa. Perovskite man- 
ganites are especially promising as they show strong 
magnetoelectronic and magnetostructural effects, includ- 
ing colossal magnetoresistance Isjl and magnetic-field- 
induced structural transitions [10|. Indeed, a first- 
principles survey of the phonon dispersions of cubic per- 
ovskite oxides, including previously reported results on 
chromites ll[ and our own investigations of chromites, 
ferrites and manganites compounds [l3] shows strong 
spin-phonon coupling in both SrMnOs and CaMnOa. 
Though both compounds have moderate a moderate Tn, 
we focus on SrMnOs as the strong oxygen octahedron 
instabilities in CaMnOa are less favorable for realization 
of the spin-phonon mechanism. 

In this paper, we report first-principles observation 
of strong spin-phonon coupling and a resulting epitaxial 
strain induced multiferroic phase in SrMnOa (SMO), a 
AFM-PE perovskite oxide. In bulk, SMO is observed 
to have G-type antiferromagnetic ordering at Tn=~233- 
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260 K in the cubic perovskite structure 
G-AFM phase, we find the epitaxial-strain-induced fer- 
roelectricity previously suggested in Ref. [isi Further, 
we find that the polar instability is strongly coupled to 
ferromagnetic ordering. Thus, the increase in polariza- 
tion produced by increasing tensile strain leads to addi- 
tional magnetic transitions, first to G-AFM and then to 
A-AFM ordering, with an increasing fraction of parallel 
nearest neighbor Mn spins, before the final transition to 
the ferroelectric FM phase. For compressive strain, the 
transition occurs directly from the G-AFM to the ferro- 
electric FM phase. While the cubic FM phase is metal- 
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TABLE L Calculated lowest phonon frequencies, in cm^^, of 
cubic SrMnOa at calculated equilibrium lattice constants with 
G-AFM and FM orderings for high symmetry g-points. 





r 


X 


R 


M 


G-AFM (ao=3.845A) 


121 


116 


84.5i 


38.1 1 


FM (ao=3.845A) 


76.2i 


116 


114? 


86.3i 


FM (ao=3.865A) 


109j 


113 


119j 


89.9i 



lie, in both cases, the system with the equilibrium polar 
distortion is insulating. Finally, we predict large mixed 
magnetic-electric-elastic responses in the vicinity of the 
phase boundaries. 

First-principles calculations were performed using 
density-functional theory within the generalized gradi- 
ent approximation GGA+C/ method with the Pcrdcw- 
Beckc-Erzcnhof paramctrization [lit implemented in 
the Vienna Ab Initio Simulation Package (VASP- 
4.6) [13, [3; selected LSDA-I-J7 calculations were per- 
formed for comparison. We use the Dudarev [l^ imple- 
mentation with on-site Coulomb interaction U=2.7 eV 
and on-site exchange interaction Jh = ^ cV to treat the 
localized d electron states in Mn. Within GGA+U, this 
choice gives agreement between the calculated {2.7 ob) 
and experimental magnetic moments (2.6/^s±0.2) 11411 
and is similar to that in a previous GGA+U study[20|. 
The projector augmented wave (PAW) potentials [2]J ex- 
plicitly include 10 valence electrons for Sr (3s23p64s2), 
13 for Mn (3p63d54s2), and 6 for oxygen (2s22p4). 

The phonon frequencies of the ideal cubic pcrovskite 
G-AFM and FM reference structures were computed us- 
ing the frozen phonon method in a x V2 x 2 supercell 
with a 6 X 6 X 4 Monkhorst-Pack (M-P) /c-point mesh 
at the r, R, X and M points of the primitive pcrovskite 
Brillouin zone; as spin-orbit coupling was not included, 
these wavevectors remained good quantum numbers for 
both FM and G-AFM orderings. 

To find the minimum-energy configuration in a given 
space group determined by freezing in one or more un- 
stable modes of the cubic reference structure, we moved 
the atoms according to the conjugate-gradient algorithm 
until the residual HcUman-Feynman forces were less than 
l.OmcV/A. Structural optimizations were performed for 
20-atom ^/2 x V2 x 2 with a 4 X 4 X 4 M-P fc-point mesh; 
for A-AFM Rj[110]-fr4 [110], a 2x2x2 superccU with a 
4x4x4 M-P fc-point mesh was used. 

To study the effects of epitaxial strain, we performed 
"strained-bulk" calculations [S^l in which calculations 
were performed for the periodic crystal with appropriate 
epitaxial constraints imposed on the in-plane lattice pa- 
rameters, with all atomic positions and the out-of-plane 
lattice constant optimized. Epitaxial strain is here de- 
fined relative to the computed lattice constant for the 
G-AFM cubic pcrovskite structure (3.845A). Ferroelec- 



tric polarizations for the relaxed structures at each strain 
were computed by the Berry-phase method [2^. Curie 
and Neel temperatures for a given strain are estimated 
from the energy differences between FM and G, G and 
A-AFM orderings assuming two exchange constants (in- 
plane and out-of-plane nearest neighbor couplings) and 
applying mean field theory As exchange couplings 
are generally overestimated in DFT, to obtain a pre- 
diction useful for quantitative comparison with experi- 
ment, we then uniformly rescaled the temperature so that 
the energy difference for the bulk cubic pcrovskite phase 
(A£'=:86meV, T^t.at ft— 3300K ) corresponds to the ex- 
perimental value of T/v=250K [H, [H. 

The calculated lattice constants of cubic G-AFM and 
FM SrMnOs are a=3.845A and 3.865A, respectively; the 
slight overestimate relative to the experimental G-AFM 
value a=3.80A 



IJI is typical of GGA calculations 
for oxides. The computed lowest phonon frequencies of 
cubic FM and G-AFM SrMnOs at the computed lattice 
constant for the cubic G-AFM structure are shown in 
Tab. [U for the cubic FM structure, the computed lattice 
constant is very similar and the effect of the difference on 
the phonon frequency, also is given in Table HI is small 
except at F. The unstable modes at the M and R points 
are the oxygen octahedron rotations. Across the Bril- 
louin zone, the FM modes are lower in frequency than 
the corresponding modes in the G-AFM structure. This 
effect is especially dramatic at F, where the lowest fre- 
quency TO mode is stable in the G-AFM structure but 
unstable in the FM structure. The corresponding eigen- 
vectors are (Sr,Mn,0|| ,0^) = (0.10,0.31,-0.52,-0.56) for 
G-AFM and (0.03,0.41,-0.42,-0.57) for FM, showing dis- 
placement of both the Sr and Mn cations relative to a 
fairly rigid oxygen octahedral network. 

From the presence of unstable modes, it is clear that 
the computed ground state of SrMnOs is not the ideal 
cubic pcrovskite structure. However, the energy gains 
and distortions for the G-AFM structure resulting from 
freezing in the dominant Mj^ and/or R^ unstable modes 
are quite small. For example, the equilibrium R^[001] 
rotation angle is 4.8°, resulting in a lA/mcm structure 
with an energy 6.2 meV/f.u. lower than the ideal cubic 
pcrovskite structure. i?3c(Rj[lll]) and /mma(Rj[110]) 
yield similar equilibrium angles and energy gains, while 
P4/m&TO(M;]"[001]) is smaller still, with 2.4° and energy 
gain 0.4 meV/f.u. This is consistent with the experimen- 
tal observation that AFM SrMnOa has a cubic structure 
but transforms to the I A/ mem, structure with a small 
fraction of ^-site substitution by the smaller cation Ca 
[l3|. For FM ordering, the energy of the equilibrium 
cubic structure is 76 mcV/f.u. higher than that of the 
G-AFM cubic structure, and all FM structures are higher 
in energy than the lowest-energy AFM structure, despite 
the lower frequencies of the M. R and F modes and larger 
relaxation energies. 

The restrictions on lattice vectors at 0% epitaxial 
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FIG. 1: (Color online) GGA+U total energies of various 
structures, per formula unit, obtained by freezing in the given 
mode(s). Calculations were performed at integer values of 
strain and interpolated. The energies of structures with FM 
ordering are shown in red, for G-AFM in blue, C-AFM in 
green and A-AFM in violet. Vertical black dotted lines at 
1.0%, 2.0%, 2.8% and 3.4% strain indicate phase bound- 
aries separating Imma{G- AFM), Jma2(G-AFM), Ima2{C- 
AFM), Ima2{A-AFM) and /ma2(FM) states. The inset is a 
schematic showing the stability of the FM-FE state at large 
strain. 



strain change these distortions and energy differences 
only shghtly; the results for the latter are included in 
Figure [TJ For example, the tetragonal relaxation of the 
FM state lowers its energy slightly, so that the AFM- 
FM splitting is 84 meV/f.u. For the R:{^[001] rotation, 
the equilibrium angle is almost the same (4.7°) and the 
energy lowering is slightly smaller (5.4 meV/f.u.). 

In Fig. [1] the epitaxial strain dependence of the to- 
tal energies of various structures and magnetic order- 
ings is shown for tensile strain from 0% to 5%. The 
energies of the G-AFM phases are seen to increase with 
increasing tensile strain. Because of their coupling to 
strain, the lowest frequency polar modes in the G-AFM 
PA/mmm and /TOma(Rj[110]) structures become un- 
stable at a critical strain of about 1%, with second- 
order phase transitions to the polar Amm2(r^[110]) 
and /TOa2(R|[110]-hr4 [110]) phases, respectively. The 
polar distortion and associated energy gain increase 
with increasing strain. This epitaxial-strain-induced 



ferroelectricity in the antiferromagnetic state is analo- 
gous to that previously found in nonmagnetic PA/mmm 
SrTiO3 0, 0], Pnma CaTiOa and antiferromagnetic 
Pnma CaMnOa [l^, with critical strains of 0.6%, 2%, 
and 2% respectively; hypothetical antiferromagnetic cu- 
bic BaMnOa [2^ is already ferroelectric at its equilibrium 
lattice constant. The dependence of the total energies on 
compressive strain, not shown, produces an analogous 
effect at -2.9% strain 

Here, we focus instead on the interplay of the strain 
and polar instability with the magnetic ordering. In con- 
trast to the AFM phases, the energy of FM structures 
initially decreases with increasing tensile strains (Fig.[T]); 
this corresponds to the slightly larger computed lattice 
constant for FM ordering in the undistorted cubic struc- 
ture. The polar instability leads to a substantial energy 
lowering which increases with increasing strain, so that at 
tensile strains above 3.4% the energy of the Ima2 FE-FM 
structure drops below those of the AFM-FE phases (see 
inset); for compressive strains, the FE-FM structure with 
space group /4cm(R4 [001]-|-r4 [001]) is favored above a 
critical strain of -2.9%. The polar distortion with strain 
also drives a metal-insulator transition above a critical 
amplitude; at all strains considered, the FE-FM phase is 
insulating, with a band gap ranging from 0.22 eV at 0% 
strain to 0.51 eV at 5 % strain. 

In addition to G-AFM ordering, we also considered 
R^[110]+r4 [110] phases with G-AFM and ^-AFM or- 
derings, in which some fraction of nearest neighbor Mn 
moments are parallel. Though higher in energy than G- 
AFM at 0% strain, they are favored by increasing tensile 
strain. As can be seen in Fig. [TJ the polar G-AFM phase, 
with 1/3 parallel-spin bonds, drops below G-AFM at 2 
%, and the A- AFM phase drops with 2/3 parallel-spin 
bonds below the C-AFM phase at 2.8 %. For compres- 
sive strain, however, the distorted G-AFM state is lower 
in energy than the A and G states up to the strain where 
the FM state becomes lower in energy, and thus there are 
no intermediate transitions. 

There are thus four phase transitions in the range of 
tensile strain considered: G- AFM-FE G-AFM-FE 
G- AFM-FE ^-AFM-FE FM-FE. The strain depen- 
dence of the electric polarization and of the estimated 
Neel (for AFM phases) and Curie (for the FM phase) 
temperatures are shown in Fig. ^ The three first-order 
magnetic transitions are of particular interest. At each of 
these, the change in magnetic order is accompanied by a 
change in the magnitude of the electric polarization. The 
magnitude of the polarization change makes typical elec- 
tric energies {P■E=10^J.C/cln'^ ■ 50 kV/cm=0.2 meV/f.u.) 
comparable to typical magnetic energies (M ■ H=Z^b • 
1T=0.2 meV/f.u.). Thus, near the A-AFM-FE^FE-FM 
phase boundary, an applied electric field can induce a 
nonzero magnetization. A substantial magnetodielectric 
coupling is also expected at that phase boundary as a 
magnetic field induced transition to the FM phase will 
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FIG. 2: (Color online) (a) Computed ferroelectric polariza- 
tion and magnetic transition temperature of SrMnOs G-AFM 
(square), C-AFM (down triangle), A-AFM (up triangle) and 
FM (circle) in the lowest energy structure at each strain value. 
The fractions of parallel nearest neighbor spins in each state 
are given. The bold black line follows the polarization of 
the lowest energy structure with changing strain. The linear 
phase diagram at the top, not on the same scale, shows the 
phases and phase boundaries as computed within LSDA+17 
for comparison. Open and solid symbols represent structures 
where R rotation is along [110] and [001] respectively, (b) 
Magnetic ordering temperatures were estimated from com- 
puted values of AE as described in the text, (c) c lattice 
parameter for the lowest energy structure at each strain value. 



also lead to a jump in the FE polarization and in the di- 
electric constant. At the AFM-FE/FM-FE phase bound- 
ary for compressive strain, in addition to an analogous 
magnetodielectric coupling, a strong strain response is 
expected due to the jump in c lattice parameter. 

As discussed above, this mechanism for epitaxial sta- 
bilization of a FM-FE phase based on spin-phonon cou- 
pling was first discussed for EuTiOa 0. In that system, 
the polar phonons for the bulk cubic phases are stable, 
with the phonon in the higher-energy FM phase about 
10 cm~^ lower in frequency than that in the G-AFM 
phase. In contrast, in SrMnOa at 0% strain both the 
FM-AFM energy difference and the FM-AFM difference 
in polar phonon frequency are much greater, the FM po- 
lar phonon being unstable. As a result, the metastable 
FM-FE structure is already present at 0% strain, though 
the electric field needed to induce this state is inaccessibly 
high. As in EuTiOs, strain acts both to increase the fer- 



roelectric instability in both phases and to lower the FM- 
AFM energy splitting. This is apparent in the lowering of 
Tn for G-AFM as the phase boundaries are approached, 
and in the increase of Tc with increasing magnitude of 
strain. The scale of the magnetic ordering temperatures 
at the phase boundaries is set by the coupling of the 
magnetic energy splitting to the discontinuous change in 
structure as the boundary is crossed; this coupling is so 
strong that the magnetic ordering temperatures remain 
above 92 K throughout the strain range considered, even 
at boundaries where the structural discontinuity is rela- 
tively small. This offers operating temperatures higher 
than that of EuTiOs and magnetically-driven fcrro- 



electrics [26|,l27|. 

In GGA+U, the overestimate of the cell volume may 
lead to a spurious enhancement of the polar instability. 
We have investigated this effect by performing calcula- 
tions for the four lowest-energy phases with LSDA-|-[/, 
where the underestimate of the volume is expected to 
lead to a comparable suppression of the polar instability. 
As can be seen from the linear phase diagram at the top 
of Fig. [21 epitaxial-strain-induced ferroelectricity in the 
G-AFM phase moves to larger strain, and the relative 
stability of the G-AFM and A-AFM phases decreases, so 
that the transition from G-AFM-FE to FM-FE occurs at 
a higher critical strain of 5.5%, which can be regarded as 
an upper bound. We thus estimate the critical s for the 
FM-FE phase to be -h4.5 ± 1% (tensile strain), and -4.9 
± 2% (compressive strain). 

It is well known that changes in the choice of the pa- 
rameter U can have a significant effect on magnetic or- 
dering energies; this was shown in particular in EuTiOs 
01 . Wc have verified that the epitaxial-strain- induced 
multiferroicity occurs for a wide range of possible choices 
of U, though the quantitative details change. For exam- 
ple, use of a slightly larger value of U reduces the energy 
difference between FM and AFM phases while leaving 
the spin-phonon coupling relatively unchanged so that 
the critical strain for the FM-FE phase shifts down from 
3.4%. 

The critical strain for observation of the FM-FE phase 
in SrMnOs is rather high, which may make experimental 
confirmation challenging. Observation of the behavior 
characteristic of lower strains, such as the strain-induced 
ferroelectricity in G-AFM and the decrease in T^v in the 
G-AFM PE and FE phases, should be taken as indicators 
of the impending transition to the FM-FE phase at higher 
strain. 

In summary, wc have presented first-principles density- 
functional calculations that reveal a large spin-phonon 
coupling in cubic SrMnOa, with ferromagnetic ordering 
producing a polar instability. Through combination of 
this coupling with the polarization-strain coupling char- 
acteristic of perovskite oxides, both tensile and compres- 
sive epitaxial strain drive the system through a series 
of phase transitions to a ferromagnetic-ferroelectric mul- 
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tifcrroic state. As the magnetic transitions are accom- 
panied by a jump in electric polarization, there is the 
possibility of electric field control of magnetic ordering; 
at the two boundaries between AFM and FM phases, 
the polarization can conversely be controlled by an ap- 
plied magnetic field. At this boundary for compressive 
strain, the jump in c lattice parameter also should yield 
a strong strain response to applied fields. Though the 
cubic FM phase is metallic, the polar distortion opens a 
gap in the electronic density of states, resulting in the 
insulating character of the FM-FE phase. This suggests 
that the search for epitaxial-strain-induced multifcrroics 
could be productively extended to include other systems 
with metallic character for the FM reference state. 
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